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METHOD AND APPARATUS FOR IMPROVING CHANNEL 
ESTIMATE BASED ON SHORT SYNCHRONIZATION CODE 

CROSS-REFERENCE TO RELATED APPLICATIONS 
[0001] This application is related to the following co-pending and commonly 

assigned patent application(s), all of which applications are incorporated by reference 
herein: 

[0002] Application Serial No. --/—,—, entitled "METHOD AND APPARATUS 

FOR REMOVING CODE ALIASES WHEN USING SHORT SYNCHRONIZATION 
CODES," filed on same date herewith, by Haitao Zhang, attorney's docket number 
020305. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

[0003] The present invention relates to systems and methods for communicating 

information, and in particular to a system and method for estimating the impulse response 
of a communication channel using short synchronization codes. 

2. Description of the Related Art 

[0004] In packet-based communication systems, spreading codes are used for 

packet detection and synchronization purposes. Correlation techniques are used to 
identify and synchronize to its timing. In many instances, the spreading code sequence 
can be in the order of 1000 chips or more. Since the receiver must correlate through all 
possible delays, this process can result in unacceptable delays. 

[0005] To ameliorate this problem, a short spreading code with good aperiodic 

autocorrelation can be used for packet detection and synchronization purposes. One 
example is the IEEE 802.11 Wireless Local Area Network (WLAN) system, which uses a 
length 11 Barker code as a spreading sequence for the preamble and the header of a 
packet. The short length of the spreading sequence makes it easy for receivers to quickly 
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detect the presence of a packet in the communication channel and to synchronize to its 
timing. 

[0006] In the case of a linear channel, for the purpose of receiver design, it is 

often desirable to estimate the impulse response of the communication channel. In the 
context of the WLAN, a multi-path linear channel is often utilized, and such 
communication channels require equalization for effective reception. Given an estimate 
of the impulse response of the communication channel, we can directly calculate 
equalizer coefficients through matrix computations, as opposed to the conventional 
adaptive algorithms. This is described in "Digital Communications," by John G. Proakis, 
4th edition, August 15, 2000, which reference is hereby incorporated by reference herein. 
This allows equalizer coefficients to be computed in a digital signal processor (DSP) 
instead of in more expensive and less adaptable dedicated hardware implementing the 
adaptation algorithms. 

[0007] Unfortunately, because the spreading code used is short (e.g. on the order 

of 11 symbols) a straightforward correlation using the spreading code will produce a 
distorted estimate. What is needed is a simple, computationally efficient technique that 
can be used to compute substantially undistorted communication channel impulse 
response estimates, even when the received signal was chipped with a short spreading 
code. The present invention satisfies that need. 

SUMMARY OF THE INVENTION 
[0008] To address the requirements described above, the present invention 

discloses a method and apparatus for estimating a communication channel impulse 
response h(t). The method comprises the steps of generating co m {t) = co{t + mNT c ) for 

m = 0,l,A ,M by correlating a received signal r(t) with a spreading sequence S i of 
length N, wherein the received signal r{t) comprises a chip sequence Cj applied to a 
communication channel characterizable by an impulse response h(t) y and wherein the 
chip sequence c } is generated from a data sequence d i spread by the spreading sequence 
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S i ; generating an estimated communication channel impulse response h M (t) as a 
combination of co m {t) and d m for m = 0,l,A,M; and filtering the first estimated 

communication channel impulse response h M it) to generate the estimated communication 
channel impulse response h(t) with a filter / selected at least in part according to the 
spreading sequence S r The apparatus comprises a correlator for generating 
co m (t) = co(t + mNT c ) for m = 0,l,A ,M by correlating a received signal r(t) with a 
spreading sequence S i of length N, wherein the received signal r(t) comprises a chip 
sequence c j applied to a communication channel characterizable by an impulse response 
h{t), and wherein the chip sequence c . is generated from a data sequence d, spread by 
the spreading sequence 5 ( ; an estimator for generating an estimated communication 
channel impulse response h M (t)as a combination of co m {t) and d m for m =0,1, A ,M ; 
and a filter /selected at least in part according to the spreading sequence 5, , the filter for 

filtering the first estimated communication channel impulse response h M (t) to generate 
the estimated communication channel impulse response h(t) . 

[0009] The foregoing permits the impulse response h(t) of the communication 

channel to be accurately estimated, even with short chip codes. Non-intuitively, in the 
case of a time-limited channel impulse response, the present invention yields an estimate 
that can be made perfect in the limit of high signal-to-noise ratio (SNR). 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0010] Referring now to the drawings in which like reference numbers represent 

corresponding parts throughout: 

[0011] FIG. 1 is a diagram of a transceiver system; 

[0012] FIG. 2 is a block diagram illustrating process steps that can be used to 

implement the present invention; 

[0013] FIG. 3 is a diagram of a transceiver system utilizing a filter / to improve 

the estimated communication channel impulse response; 
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[0014] FIG. 4 is a diagram showing the response of the filter; 

[0015] FIG. 5 is a flowchart describing exemplary processing steps that can be 

used to improve the reconstruction of the value of the communication channel impulse 

response using super codes imposed on the portion of the data sequence; 

[0016] FIG. 6 is a diagram of a transceiver system utilizing super code to transmit 

sequences; 

[0017] FIG. 7 is a diagram presenting a correlator output using 11 symbol long 

Barker code; 

[0018] FIG. 8 is a diagram presenting a correlator output using Walsh codes as an 

input super code; 

[0019] FIG. 9 is a diagram presenting a correlator output after postprocessing with 

a filter /as described in FIG. 2 and FIG. 3; 

[0020] FIG. 10 is a diagram presenting a more detailed view of the main lobe 

peak, showing the estimate of the communication channel impulse response in the actual 
communications channel impulse response; and 

[0021] FIG. 1 1 is a diagram presenting one embodiment of a processor that can be 

used to practice the present invention. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
[0022] In the following description, reference is made to the accompanying 

drawings which form a part hereof, and which is shown, by way of illustration, several 
embodiments of the present invention. It is understood that other embodiments may be 
utilized and structural changes may be made without departing from the scope of the 
present invention. 

System Model 

[0023] FIG. 1 is a diagram of a transceiver system 100. Using signal spreader 103, 

a random data symbol sequence d t 102, comprising a series of data packets 128 (each of 
which may include a preamble 124 used by the receiver for identification purposes, as 
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well as a data payload 126) is spread by a sequence S g 104 of length N: {S n ,0<n<N-l\ 
and having a chip period. The sequence S. 104 is known to the receiver 112 apriorL The 
spread chip sequence c y 106 is therefore: 

Cj = c iN+n =d i .S n ,0<n<N-l Eq. (1) 

[0024] This spread chip sequence c . 106 is transmitted through a linear 

transmission channel 108 having a combined channel impulse response h(t). The 
transmitted signal is received by a receiver 112. The received waveform r(t) 1 14 is: 

K0= Zs*/i(r- 7 T c )-hn(0 Eq.(2) 
where n{t) 121 is an additive noise component. 

[0025] This formulation does not explicitly impose a causality requirement on 

h(t) 108. If explicit causality is desired, this can be accomplished by setting 
h(t) = 0j < 0. For simplicity purposes, all the data and code sequences in the following 
discussion are assumed to be real, though the channel impulse response h(t) 108 and the 
additive noise component n(t) 121 could be complex in their baseband representations. 
Complex sequences could be easily accommodated if needed, but they are not common 
for synchronization purposes. 

[0026] The receiver 112 receives the transmitted signal, and correlates the 

received signal r(t) 114 with the known spreading sequence 5, 104 to identify the data as 
intended to be received by the receiver 1 12. Once the received signal r(t) 1 14 is received, 
the preamble can be examined to determine the address of the data and whether further 
processing is necessary. 

[0027] Such systems also use the received signal to estimate the input response of 

the communication channel 108. This information is used to improve later detection and 
reception of signals from the transmitter 110. In circumstances where the spreading 
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sequences, 104 is relatively short, the data packet 128 must be detected quickly, and 
there is less data available to estimate the response of the communication channel 108. 

Conventional Detection and Synchronization 
[0028] For detection and synchronization purposes, the search for the spreading 

code is conventionally performed by correlating the received signal r(/)114 with the 
spreading sequence. This is accomplished by the correlator 116. Although this 
correlation is typically done after sampling in the time domain, for notational simplicity, 
we do not perform the time domain discretization. The correlator 116 output co(t) 118 is 
given by: 

AM 



co(t) = ]>>(* + (#- lK - iT c ) • S„_,,, Eq. (3) 

= £r(f+iT e ).5, Eq.(4) 

1=0 

= 1 1>, •h{t-(j-i)T c )»S i +n{t) Eq. (5) 

= f i D{l).h(t-lT c )+n{t) Eq.(7) 

where D(l) is the correlation between the chip sequence and the spreading sequence and 
we will refer to it as the chip correlation. 

[0029] For notational simplicity, we have introduced a (negative) group delay 

(IT C ) in calculating the correlator output 118. The correlator 116 output is given by the 
convolution of the chip correlation D(l) with the sampled communication channel 
impulse response hit - IT C ) plus a noise component n (t) . Upon further examination: 

£>(0 = Zc (+ ,*5, Eq.(8) 
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= Z d n, * S ^ • S, , + £ rf m+1 • S n+i _„ • 5,- Eq. (9) 

j=0 i=N-n 

= d m .A(n) + d m+l .A{N-n), Eq. (10) 

l = mN + n,0<n<N Eq. (11) 

where A(n) is a two-sided aperiodic autocorrelation of the spreading sequence defined as: 

A(n)=A(-n)= . S i+n ,0<n< N 

1^ Eq. (12) 

A(n) = 0,|rc| > N 

A(n) is a property of the code sequence that is known by the correlator 116 

apriori. 

[0030] For detection and synchronization purposes, the spreading sequence S i 

104 is designed to have minimum values of A(k) when k ^ 0 . However, for small (e.g. 
on the order of 10) values of N (short spreading codes), even the smallest side lobe 
magnitude is not negligible compared to the in-phase autocorrelation. 
[0031] Barker sequences, when they exist, give the best aperiodic autocorrelation. 

For an 11 chip Barker sequence, 5. = 1,-1,1,1,-1,1,1,1,-1,-1,-1 , the autocorrelation 

becomes A(i) = 11,0,-1,0,-1,0,-1,0,-1,0,-1 for 0 < £ < 1 1 . Note that even for Barker 
codes, because the spreading sequence S i 104 is of limited length, the 
autocorrelation A(i) includes significant side lobes. 
[0032] The correlator 1 16 output 118 can be rewritten as: 

co{t)= XX^OW + i).h{t-{jN + i)T c )+n(0 Eq.(13) 

j=-oo ,=0 

= Z •Mfhdj.t •A(N-i))*h{t-(jN + i)T c ) + n{t) Eq. (14) 

= Z IX • # h (* ~ U N + i)T c )+A{N-i). h{t - {{j -l)N + i)T c )) + n(t) 

Eq. (15) 
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= X ^dj*Aii)*h{t-jNT c -iT c )+n{t) Eq. (16) 

= idj>h{t--jNT c )+n{t) Eq.(17) 

where the following is defined as the convolution of the spreading sequence aperiodic 
autocorrelation A(i) and the sampled channel impulse response h(t - iT c ) as follows: 

*(0= I>(i>M'-*T e ) Eq.(18) 

[0033] This is an estimate of the combined communication channel 108 impulse 

response h(t) at the output of the code correlator 116. 

The above equations can be more succinctly written using a convolutional 
notation. Defining a convolution of two infinite sequences A i and B i as 

C = A®B&C(i)=^ A (j) • B(i - j\ V/ Eq. (19) 

y 

[0034] By defining an operator O that converts any sequence O to a time domain 

function using the Dirac delta function: 

B T {t) = J^B(i)*S(t-iT) Eq.(20) 

i 

we can also define the convolution of a function with a sequence using a normal 
convolution of two functions: 

C{t) = A{t)®B T <=>C{t) = ^A{t-jT)B{j) Eq.(21) 

j 

[0035] Using the above notations and further, by adopting the following 

definitions: 

u(/jV) = rf. (data) Eq.(22A) 
u(iN + n) = 0fi<n<N Eq. (22B) 

f S n ,0<n<N 

n) = \ (a time-limited chipping sequence) Eq. (22C) 

^0, otherwise 

the foregoing equations (1), (2), (3), (6), (12), (18), (16), (17) can be rewritten as: 
c = u®S Eq.(l') 
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r = h®c' +n o 


Eq. (2') 


co = r® Sl r 


Eq. (3') 


= h®c T '<8> Sl r + n, where 5_ («) = 5(- n) 


Eq. (6') 


A = S®S_ 


Eq.(12') 


h = h®A T < 


Eq.(18') 


co = /i ® m c +n 


Eq.(16') 


= h®d m < +n 


Eq. (17') 



Determining a Communication Channel Impulse Response Estimate 
[0036] For simplicity of notation, in the remaining discussion, we assume that 

data symbols are binary. The results however, can be generally applied to non-binary 
data. 

[0037] Because the correlator 116 has access to the same code sequence S i 104 

that was used to generate the spread chip sequence Cj 106 before transmission, the 

correlator 116 can correlate the received signal r(t) 114 with the code sequence S,. 104. 

However, aliasing can occur with short code sequences 5, 104, because time delays may 

cause the correlator 116 to correlate different portions of adjacent code sequences. 
Conventionally, these aliasing effects are reduced by integrating or summing over 
multiple (e.g. M) code periods, as discussed below. 

[0038] As described in Eqs. (13)-(17), based on the correlator 116 output 118 we 

can form an estimate of the channel impulse response over one code period T c : 

h x (t) = d 0 »co(t) 

= h{t) + X do • d } • h(t - jNT c ) + d 0 • n(t) Eq * (23) 

where d 0 is a value of the data at time t = 0 . 
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[0039] This is a rough approximation to h(t), corrupted by aliased copies of h(t) 

spaced at multiples of NT C away from the desired copy. These aliasing and the additive 
noise terms can be reduced through further summation over M code periods: 

K (t) = XT IX * c °(' + mNT c ) E( i- C 24 ) 



M m=0 

A/-1 



= ) + tV Z Z • <*y # + ( w ~ ^ ) + »* (0 Eq- ( 25 ) 

= *(0 + ^rZZ rf « + % (0 Eq. (26) 

[0040] The foregoing indicates that through output 122 of estimator 120, by 

removing the data modulation through correlation with the data sequence, we obtain an 

estimate h of the channel impulse response plus terms defined by the autocorrelation of 
the data sequence, which vanish when summed over infinite terms. 
[0041] If D M (/) is defined as: 

D *( l ) = ^Ild m .d Um Eq,(27) 

then 

h M =h®D^+n' M Eq.(28) 
wherein h M is an estimate of the communication channel impulse response h(t) . 

When the data sequence d i 102 is random, white and independent of the additive noise 

n(i) 121, and in the limit of M -» °°: 

A. (0 = =° 

r r Eqs. (29) 

h„=h 

[0042] Therefore, in the limit of infinite summation (as M approaches infinity), 

we obtain an estimate that is equal to the true channel impulse response h(t) convolved 
with the aperiodic autocorrelation of the spreading sequence S t 104. 

[0043] As the foregoing demonstrates, we can not obtain the true channel impulse 

response h(t) with simple integration. The best we have is smeared by the autocorrelation 
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of the spreading sequence 5. 104. In cases where the spreading sequence S { 104 is long, 
the autocorrelation approaches a delta function, and the side lobes disappear. However, 
when the spreading sequence 5, 104 is short, the sidelobes of the autocorrelation are not 

negligible and will cause significant distortion to the estimate of the communication 
channel impulse response h(t) . 

Improved Channel Estimates for Short Spreading Sequences 
[0044] As is demonstrated below, the present invention improves the 

communication channel impulse response estimate by filtering the first estimated 
communication channel impulse response h M (t) to generate the estimated communication 
channel impulse response /i(r)with a filter /selected at least in part according to the 
spreading sequence 5 ( . In particular, when the time span of the communication channel 

108 is limited, a zero-forcing deconvolution can be used to improve the estimate. 

[0045] FIG. 2 is a block diagram illustrating process steps that can be used to 

implement the present invention. 

[0046] FIG. 3 is a diagram of a transceiver system 300 utilizing the filter / 

described above to filter the first estimated communication channel impulse response 
h M (t) to generate an improved estimate suitable for short spreading sequences S i 104. 
[0047] Referring to FIG. 2 and FIG. 3, blocks 202 through 208 recite steps that 

are used to generate co m (t) 118. A spread chip sequence Cj 106 is generated from a data 

symbol sequence d { 102 and a spreading sequence S i 104 of length N, as shown in block 
202. The chip sequence Cj 106 is transmitted via a communication channel 108 as 

shown in block 204, and received as shown in block 206. The communication channel 
includes the transmitter 110 and the receiver 112. The received signal r(t) 114 is then 
correlated with the spreading sequence 5, 104, by the correlator 116 to generate co m {t) 
as shown in block 208. 
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[0048] In block 210, an estimated communications channel impulse response h M (t) is 

generated by the estimator 120 as a combination of co m {t) and d m for m = 0,l,K ,M 

This can be accomplished, for example using the relationship described in Eq. (24) above. 

[0049] Finally, in block 212, the first estimated communication channel response h M (t) 

is filtered with a filter /selected at least in part according to the spreading sequence S ( . 

104. In one embodiment, the filter is a finite impulse response (FIR) filter /302 designed 
with the following constraints: 

A f ^A®f Eq.(29) 
A f (0) = 1, A f (n) = 0,0 < \n\ < L Eq. (30) 

wherein A ® f is the convolution of the autocorrelation of the spreading sequence S { 
104 and the filter, and A f is the autocorrelation of the spreading sequence S f 104 after 
filtering. 

[0050] FIG. 4 is a diagram showing the response of the filter/ 302 described in 

Eqs. (29) and (30). 

[0051] When the estimate of the communication channel impulse response is 

filtered with this filter, we obtain: 
h f =h®f T < 

T T 

= fc®A 7; ®/ 7; Eq. (31) 

= h® Af 

[0052] Using this technique, the effects of the side lobes (aliased versions of the 

autocorrelation of the spreading sequence S i 104) are eliminated between L and -L. The 

side lobes are not completely removed (since the filter passes components greater than L 
and less than -L) but the result near the origin (n = 0) is of primary interest, and the effect 
of the side lobes can be significantly reduced in this region. 

[0053] If the time span (duration of the impulse response) of the communications 

channel is less than LT C , i.e. 
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3'i <t 2^2~ t i <LT c yt<t } \Jt>t 2 :/z(r) = 0 Eq. (32) 

(that is, there exists a time t 2 greater than /, defining a time interval t 2 - t l less than LT C , 
and for all time outside of the interval t 2 -t x , h(t) is close to zero), 

[0054] Then, the filtered estimate h f (or, in the earlier notation, h f (t)) is 

composed of an exact copy of h(h(t)), plus some aliased versions of it in non- 
overlapping locations. So in this case h is resolvable from h f . 

[0055] Such a filter with length 2L + lcan be designed with the simple zero- 

forcing criteria: 

L 

2 Mn - i) • f(i) = A f (n) < n < L Eq. (33) 

i=-L 

wherein / (i) is the impulse response of the filter/ 302 such that A f (rc) is a convolution of 
A(n)and f(J),A f (n)=l for n = 0 and A f (n) = 0 for0<|n|<L, and 

N-l-n 

A(ri) = A(-n) = ^ S t • S i+n ,0 < n < N , and wherein Wis a length of the chip sequence 

S { 104. L can be chosen such that the product LT C (the chip period T c is known) is 
approximately equal to the time span (e.g. the approximate duration of the impulse 
response) of the channel 108. 

[0056] Note that the value A(n-i) is well defined ... it is a property of the 

spreading sequence 5. 104, which is known apriori. 

[0057] As usual, the matrix structure of the linear equations is Toeplitz. By the 

design requirement of the spreading sequence S t , the matrix should be well conditioned. 

The filter coefficients can be computed offline given the spreading sequence and desired 
window width L . 

[0058] While the foregoing has been described in respect to non recursive filters, 

other filters, such as recursive filters may also be used. A recursive filter, for example, 
may provide perfect filtering of the sidelobes, but the result may not be the quell 
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conditioned matrix, hence the solution may be more difficult to determine. In fact, any 
filter of length 2L + 1 can be defined. 



Super Coded Transmit Sequences 

[0059] It has been shown that given h and with filtering, it is possible to recover 

the true channel impulse response for a time limited channel. However, in the foregoing 

discussion, h was obtained through integration over multiple spreading sequence periods. 
The number of periods we need to integrate over can be large especially if 2L > N , since 
we rely on the autocorrelation of the data to suppress the aliased copies of h . 
[0060] In one embodiment of the present invention, supercodes, such as Walsh- 

like supercodes, are used to drastically reduce the amount of the integration required. 
This technique is especially useful in systems having sufficient a signal-to-noise ratio 
(SNR). 

[0061] Consider a pair of length 2 Walsh codes w 0 = {+ l,+l} and w l = {+ l,-l}. 

These codes can be used to form a data sequence: 
...+,+,+,---... 

[0062] Any length 2-symbol length segment from this sequence can be described 

as either w 0 or - w 0 , except for a single w x in the center. If this sequence is now 

correlated with Wj , the resulting correlation will be characterized by a single peak in the 
center and zeros elsewhere (except near the boundaries). Negatives of the two codes may 
be taken (e.g. w 0 = {-l,-l} and w, = {-l,+l} and/or their roles may be swapped (e.g. 

w, ={+l,+l} and w 0 ={+l,-l} with the same result. The three additional patterns thus 

obtained and their correlator patterns are listed below: 
..-,--,-,+,+,+,+... 
...-,+,-,+,+,-,+,-... +,+ 
-j 1 1 h — 



Express Mail No. EL 977101 138 US PATENT 
020306 

-15- 

[0063] Since the following results are equivalent for all of the above patterns 

when the additive noise is uncorrelated at sampling points, we limit our discussion to the 
first data sequence (i.e. ...+,+,+,-,-,-.. ). In this case, 

d x =+l,V/, </<0 Eq.(34) 

d { =-l,V/ 2 >i>0 Eq.(35) 

h 2 (t) = - (rf 0 • co(t ) + d, • co(t + NT C )) Eq. (36) 
2 

= *(f) + Efa -d M ).h(t-lNT c ) + n&) Eq.(37) 

= *(/) + 2( rf <-^ + ,)-^-^ c )+/i 2 , (0 Eq.(38) 

lZtlUl>l 2 



[0064] If the condition that - ^TV > (2// + L) I / 2 W > (2N + L) can be satisfied, h 

can be reconstructed free of aliasing interference, and by deconvolution (aforementioned 
filtering technique), h can be reconstructed as well. 

[0065] From the foregoing, it can be determined that a small supercode imposed 

on a portion of the data sequence can provide an alias free estimate of the communication 
channel impulse response when the channel response is time-limited. The only source of 
distortion from this estimate comes from the additive noise, which can be suppressed by 
the spreading gain times a factor of 2 (to account for the supercode). When the noise is 
low, such an approach is preferable over long integrations. 

[0066] For moderate values of L, such code sequences can be easily embedded 

within a longer preamble to packet data, probably with multiple copies, without adversely 
affecting the spectrum properties of the transmission. In addition, when the signal to 
noise ratio (SNR) is low, traditional integration as outlined in the first half of this section 
can still be carried out on such a preamble to obtain a higher processing gain against the 
additive noise. 

[0067] FIG. 5 is a flowchart describing exemplary processing steps that can be 

used to improve the reconstruction of the value of the communication channel impulse 
response by using supercode imposed on a portion of the data sequence. 
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[0068] FIG. 6 is a diagram of a transceiver system 600 utilizing super coded 

transmit sequences to generate an improved communication channel impulse response 
estimate suitable for short spreading sequences S, 104. 

[0069] In block 502, a data sequence d { 102 is generated. The data sequenced. 

102 includes one or more data packets 128, each data packet having a preamble 124 
including a constrained portion Cd i 602. The preamble 124, can be, for example, in the 
form of a pseudorandom code. 

[0070] The constrained portioned, 602 is associated with at least two codes, w Q 

and Wj. The codes w 0 and w l are selected such that the correlation A code (k) of the 
constrained portion Cd { 602 and at least one of the codes w Q and w x , is characterized by 
a maximum value at k = 0 , and they value less than the maximum value at k * 0 . 
[0071] Ideally, the correlation A code (k) of the constrained portion Cd x 602 is an 

impulse, with A code (k) equal to one at k = 0 , and equal at all other values for k. 
However, because such correlation characteristics are typically not realizable, codes w 0 
and w, can be chosen to approximate this ideal. For example, codes w 0 and w l can be 
chosen such that the correlation A code {k) of the constrained portion Cd i 602 and at least 
one of the codes w 0 and w,, is such that A code (k) = l at £=0and A code (k) ~ Ofor 
substantially all k * 0 . Or, codes w 0 and w { can be chosen such that the correlation 
A code (k) of the constrained portion Cd { 602 and at least one of the codes h> 0 and w p is 
such that A code (k) = 0 for 0 < \k\ < J , wherein J is selected to minimize the correlation of 
the constrained portion Cd t with the one of the codes w 0 , w x for substantially all k * 0 . 
[0072] In one embodiment, the constrained portion Cd { 602 comprises the pair of 

length two Walsh codes in the first sequence described above. Other embodiments are 
envisioned in which the codes are of another length (other than length two), or are codes 
other than a Walsh code. 
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[0073] In block 504, a chip sequence c- 106 is generated. The chip sequence c ■ 

106 is generated by applying a spreading sequence S { 104 of length N and having a chip 
period T c to the data sequence d i 102. 

[0074] This spread chip sequence c. 106 is transmitted through a linear 

transmission channel 108 having a combined channel impulse response h(t). The 
transmitted signal is received by a receiver 112, 

[0075] In block 506, the receiver 112 receives the transmitted signal, and 

correlates the received signal r(0114 with the known spreading sequence S, 104 to 
identify the data as intended to be received by the receiver 112. This is accomplished by 
generating co m (t) = co(t + mNT c ) for m = 0,l,A,M, using techniques analogous to 
those which were described above. 

[0076] In block 508, an estimated communication channel impulse 

response h M (t) is generated as a combination of the correlation co m (t) and the data 
sequence d m for m = 0,1, A , M . 

[0077] In one embodiment, the codes w 0 and w>, are two symbol-long Walsh 

1 M ~ l 

codes, andh M (t) computed as — y\d m •co(t + mNT c ) , with M = 2. In this case, h M (t) 

M m=0 

equals h 2 (t ) = ~ (d 0 • co(t ) + d x • co(t + NT C )) . 

[0078] Hence, where the data has been constrained with a symbol such as a Walsh 

super code, an improved estimate of the communications channel impulse response can 
be obtained by taking two consecutive values of the correlation of the received data and 
the spreading sequence and multiplying each result by the data sequence. In the example 
of Walsh codes w> 0 ={-l,-l} and w x ={-l,+l} applied to the sequence... +,+,+,-,-,-... , 

and m>, applied at the receiver, the result is that one of the values of co(t) is multiplied by 
a one, and the other is multiplied by a minus one. Hence, the output will produce 
essentially no response until the transition between the two Walsh codes occurs, at which 
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time a clean, alias-free copy of the communications channel impulse response will be 
produced. 

[0079] A length 2 supercode for improved alias suppression has been described. 

When the SNR is low and longer integration period is desirable, it would appear attractive 
to generalize the code to longer lengths. Counterintuitively, this is not possible. This 
result is shown below, by presenting a definition of such codes and showing that no such 
codes with length larger than 2 exist for binary data sequences. 

[0080] An infinite sequence A forms an impulsive correlation pair with a length L 

finite sequence B if A satisfies the following equations: 
A(i) = B{i ),V0 < i < L 

^A(i + n)^B(i) = 0yn±0 

i=0 

[0081] By contradiction, it can be shown that for binary sequences, such a pair 

does not exist for L > 2 . Supposing such sequences exist, it is apparent that L must be 
even. Considering two such cases (L = 4k and L = 4k+2) 
[0082] In the first case, L = 4k, consider the first constraint: 



£A(i-l).fi(i) = 0, Eq.(39a) 

i=0 

A(- 1) • B{0) + ^ B(i - 1) • B{i) = 0 Eq. (39b) 



[0083] Since there are 4k summands in the equation taking values from {+1,-1} 

half of them or 2k terms must be positive, and the other half negative. The product of all 
the summands must therefore be 1. 



A(-i).B(o).n(*(*-i)»*(0)=i 

A(-l)«B(L-l) = l Eqs. (40) 

M.-i) = b(l-i) 

[0084] Similar arguments can be used to show that: 

A(0 = B(L + i\-L <i<0 Eq. (41) 

[0085] But this implies that: 
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£ A(i - L) . B{i) = A(- L) • 5(0)+ g A(/ - L) . fl(0 

= A(- L) • B(0) + ]T ) • «(0 E q s - ( 42 ) 

1=1 

= A(- Z>£(0)+L-1 
>0 

which contradicts the assumption that the cross-correlation is zero everywhere except at 
the origin. Hence, by contradiction, we have shown that for binary sequences, such a pair 
does not exist for L > 2 . 

[0086] A similar argument can be applied for the second case, L=4k+2, except 

that the product of all the summands in each equations must be -1, since now we must 
have 2k+l negative terms. This leads to: 

A(i) = (- 1) 1 B(L + i)-L <i<0 Eq. (43) 

[0087] When k>0, 

L-l L-\ 

£ A{i - 2) • B(i) = A(- 2) • B{0) + A{- 1) • + £ A{i - 2) • B(i) 

i=0 /=2 

= 5(L - 2) • 5(0) - B{L - 1) • + Yj B (* " 2 ) # 5 (0 

i=2 

= 0 

t-1 t-3 



i=0 i=0 
L-3 



£ A(i - L + 2) • B(i ) = 2 A(i - L + 2) • B(i ) + A(o) • B(L - 2) + A(l) • B(L - 1) 

L-3 

= X (- 1 )' B ( f + 2 ) • B (0 + #(°) • b{l - 2) + b(i) • b{l - 1) 

i=0 

L-l 

= B{L - 2) • B(0) + B(L - 1) • B(l) + £ (- l)' B(/ - 2) • B(i ) 



i=2 

= 0 



Eqs. (44) 



[0088] Summing the two equations together we have: 

2k 

B{L - 2) • B(6) + X B{2i - 2) • S(2i) = 0 Eq.(45) 



i=i 
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[0089] However, this result is clearly impossible since there are an odd number of 

terms on the left. By contradiction it is therefore shown that it is impossible to satisfy the 
constraints when L>2 for binary sequences. 



Noise Effects 

[0090] The foregoing has demonstrated that distortions due to this spreading 

sequence design can be removed from the estimate of the communications channel 
impulse response. Attention is now turned to the remaining distortion caused by the 
additive noise. n{t) 121. Assuming that the noise source is white and stationary and is 
filtered by a receiver filter for bandwidth matching, its distortion measure can be defined 
as follows: 



A = E 



= n 



L 



4ti^ <fi (' + mNT < - /7 > /(o 

= 17 * + + iT c - IT C )• 5(«). /(/) 

l M-l 

= -^YLd m *n{t + mm c + jT c ) • R fi {j) 



m=0 j 



where 



Eq, (46) 
Eq. (47) 



Eqs. (48) 



Eq. (49) 



[0091] The ensemble expectation of Eq. (46) can be taken over n(t) y whose 

autocorrelation can be determined by the front end receive filter, and is assumed to be 
known). 
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Rjr) = *k(t)'n 0 (t + T)] 

A = T7T I ZEE*- • Kjs 0> R m ((-' - m)NT c + (/ - /fc ). ^ (/)• rf„ 

= EZ^0>^((/-yX)*^(/) 



Eqs. (50) 

00 = ~ • ^ 0 + mN ) E( i- < 51 > 



[0092] When the noise n(r) is white, we have: 

RjkT c ) = 0,Vk*0 
A = /U0)l|*£(7f 



Examples 

[0093] FIG. 7 through FIG. 10 are diagrams illustrating the performance 

improvements achieved by application of the present invention. These illustrated 
examples and places them, whereby a length 11 Barker code is used as the spreading 
sequence 5. 104. FIGs. 7-10, normalized magnitudes as a function of chip timing. No 

adjustments were made for group delays introduced by correlation, filtering and 
windowing, therefore time coordinates should be treated in the relative sense. FIGs. 7-10 
also do not include the effects of additive noise. 

[0094] FIG. 7 is a diagram presenting a correlator 116 output using a length 11 

Barker code and conventional communication channel impulse response estimation 
techniques. The correlator 116 output shows a main lobe peak 702, and multiple spurious 
peaks 704. These spurious peaks 704 (which are 1 1 chips, or NT C seconds, apart due to 

the length 11 Barker code) are due to the repeated transmission of the short code S, 104, 
which are "aliased" back upon each other. If the length of the periodic spreading 
sequences, 104 were longer, there would be fewer spurious peaks 704, and the peaks 704 
would not overlap the main lobe peak 702 as much as is shown in FIG. 7. 
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[0095] FIG. 8 is a diagram presenting a correlator 1 16 output using the Walsh 

codes in conjunction with the supercode technique described in FIG. 5. To generate this 
plot, the input data was constrained with two symbol-long Walsh codes w 0 and w l , and 

the output was processed by summing two successive outputs of the correlator 116 as 
shown in Eq. (36). For the 1 1 chips on either side of the main lobe peak 702, there is zero 
correlation, and many of the spurious correlator peaks 704 that were apparent in FIG. 7 
are no longer evident. Note, however that since only six bits of the data sequence are 
constrained ...+,+,+,-,-,-..., some aliased versions of the main lobe peak 704 (labeled 
802) are present (33 chips from the main lobe peak 702). However, since these aliased 
versions 802 are widely separated from the main lobe peak 702, an accurate estimate of 
the communications channel impulse response can be obtained. Note that a similar result 
can also be achieved without constraining the input sequence with the super code, but 
this would require integration over large number of symbols (e.g. M in Eq. (26) would be 
large). Note also that the main lobe peak 704 still includes minor peaks because the 

estimator 120 produces h 9 which is a smeared version of h. These undesirable 
components 804, caused by the autocorrelation of the spreading sequence 104, cannot be 
removed by constraining the data sequence. Instead, these undesirable components 804 
can be removed by filtering as described with respect to FIG. 9 below. 
[0096] FIG. 9 is a diagram presenting a correlator 116 output shown in FIG. 8 

after postprocessing with a filter /as described in FIGs. 2 and 3. Note that the sidelobes 
802 shown in FIG. 8, have been pushed away from the main lobe peak 702, and some of 
the undesirable components 804 of the main lobe peak 702 have been filtered. Also note 
that the data indexing (the chips shown as the time axis) of FIG. 9 has changed relative to 
the data indexing of FIG. 8. As described above, this difference is an artifact of the 
software used to plot FIG. 7- FIG. 1 1 and is not associated with the applicant's invention. 
[0097] FIG. 10 is a diagram presenting a more detailed view of the main lobe 

peak 702, showing the estimate of the communication channel impulse response 
(indicated by the asterisks) and the actual communication channel impulse response. 
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Note that the estimated communication channel impulse response follows that of the 
actual response very closely. 

Hardware Environment 
[0098] FIG. 11 and is a diagram illustrating an exemplary processor system 1102 

that could be used in the implementation of selected elements of the present invention 
(including, for example, portions of the transmitter 110, the receiver 112, the correlator 
1 16, the estimator 120, or the filter 302). 

[0099] The processor system 1102 comprises a processor 1104 and a memory 

1106, such as random access memory (RAM). Generally, the processor system 1102 
operates under control of an operating system 1108 stored in the memory 1106. Under 
control of the operating system 1108, the processor system 1102 accepts input data and 
commands and provides output data. Typically, the instructions for performing such 
operations are also embodied in an application program 1110, which is also stored in the 
memory 1106. The processor system 1102 may be embodied in a microprocessor, a 
desktop computer, or any similar processing device. 

[0100] Instructions implementing the operating system 1108, the application 

program 1110, and the compiler 1112 may be tangibly embodied in a computer-readable 
medium, e.g., data storage device 1124, which could include one or more fixed or 
removable data storage devices, such as a zip drive, floppy disc drive, hard drive, CD- 
ROM drive, tape drive, etc. Further, the operating system 1108 and the application 
program 1110 are comprised of instructions which, when read and executed by the 
computer 1102, causes the computer 1102 to perform the steps necessary to implement 
and/or use the present invention. Application program 1110 and/or operating instructions 
may also be tangibly embodied in memory 1106 and/or data communications devices 
1130, thereby making an application program product or article of manufacture according 
to the invention. As such, the terms "article of manufacture," "program storage device" 
and "computer program product" as used herein are intended to encompass a computer 
program accessible from any computer readable device or media. 



Express Mail No. EL 977101 138 US PATENT 
020306 

-24- 

[0101] Those skilled in the art will recognize many modifications may be made to 

this configuration without departing from the scope of the present invention. For 
example, those skilled in the art will recognize that any combination of the above 
components, or any number of different components, peripherals, and other devices, may 
be used with the present invention. For example, an application-specific integrated circuit 
(ASIC) or a Field-Programmable Gate Array (FPGA) can be used to implement selected 
functions, including the correlator 116, and filtering functions can be performed by a 
general-purpose processor, as described above. 
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Conclusion 

[0102] This concludes the description of the preferred embodiments of the present 

invention. The foregoing description of the preferred embodiment of the invention has 
been presented for the purposes of illustration and description. It is not intended to be 
exhaustive or to limit the invention to the precise form disclosed. Many modifications 
and variations are possible in light of the above teaching. It is intended that the scope of 
the invention be limited not by this detailed description, but rather by the claims appended 
hereto. The above specification, examples and data provide a complete description of the 
manufacture and use of the composition of the invention. Since many embodiments of 
the invention can be made without departing from the spirit and scope of the invention, 
the invention resides in the claims hereinafter appended. 



